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Abstract. I review the evidence for the importance of feedback from massive 
stars at small and large scales. The feedback mechanisms include accretion 
luminosity, ionizing radiation, collimated outflows, and stellar winds. The good 
news is that feedback doesn't entirely prevent the formation of massive stars, 
while the bad news is that we don't know what does limit their masses. Feedback 
from massive stars also influences their surroundings. I argue that this does not 
produce a triggering efficiency above unity, nor does it prevent lots of prompt 
star formation in GMCs, though it may preserve massive remnants of the clouds 
for many dynamical times. 



1. Small Scale Feedback and the Upper Mass Limit 



Massive stars influence their surroundings through a wide variety of feed back 
mechanisms. Their ionizing radiation (e.g., Vacca. Garmanv. ShuU| [l996) and 
line-driven winds ([Castor. Abbott. &: Klein 19751 ) have been known for many 
years, while only more recently have collimat ed jets and wide-angle ou tflows 
been observed from young, massive stars (e.g.. iBeuther &: Shepherdll2005l ). 

The initial mas s function for s tars above a solar mass or so follows a power 
law with slo pe -2.35 ( Salpeterlll955l ) up to an upper cutoff mass of m up ~ 100 M Q 
( Figer [20051 ) . The two most favored ideas currently for what determines the slope 
of the IMF are turbulent f ragmentation of the pa rent cloud ( Padoan et al.l l2007) 
or competitive accretion ( Bonnell fe Bate! I2006I ). What determines the value 



of the cutoff mass m up remains unknown. Limitation of accretion by increas- 
ingly strong feed back in increasingly massive stars remains a viable candidate 
mechanism ( Zinnecker &: Yorke 20071 ). 

In fact, the limitation of accretion by radiative feedback had long been 
a severe problem in theories of massive star formation, because in models of 
spherically symmetric accretion, radiation pressure on stan dard interstellar dust 



preve nts accretion on to stars with mass M > 30 M ([Wolfire Sz Cassinelli 



119871 1. However, recent models have called this conclusion into question because 
of better treatment of dust properties, non-spher ical accretion, and t he high 
accretion rates expected in high pressure cores (e.g-. lMcKee fc Tanll2003l ) such as 
those that have been identified by t he Midcourse Space Experiment and Spitzer 
Space Telescope ( Shirley et al.1 [20031 ). 
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Yorke &z Sonnhalterl ( 20021 ) performed two-dimensional, nested-grid, numer- 



ical models of the collapse of slowly rotating, unmagnetized, molecular cores, 
using a wavelength-dependent treatment of dust grain properties for a mixture 
of small amorphous carbon grains, astrophysical silicates, and, in colder regions, 
water and ammonia ice coatings. They found that a combination of decreased 
opacity, disk accretion, and focusing of ra diative losses through the pol ar cavity 
(the "flashlight effect" first discussed by Yorke fc Bodenheimerl Il999l ) allowed 
formation of 35-45 M stars. Krumholz. Mckee &: Kleinl ( 2005a! ) used computa- 
tions of radiative transfer to demonstrate that the flashlight effect can be further 
enhanced if a collimated gas outflow clears a cavity in the envelope. 

The efficacy of r adiation in preventing accre ti on is further reduced by two 
other effects. First, Krumholz. Klein &; McKee ( 20051 ) showed that in three 
dimensions, the radiation-driven bubbles that form when the radiation pressure 
begins to repel g subject to Ray leigh- Taylor instability. Essentially, the 

radiation acts as a light fluid that is trying to support the heavy gas. As a 
result, the bubbles tend to collapse, allowing add itional accre tion. Second, in 
a magnetized core, the ph oton bubble instability (lAronsI 119921 ) can also act to 
allow radiation to escape (jTurner. Quataert. k, Yorkd 120071 ) . This occurs when 
compressive MHD waves are amplified by radiation pressure, forming low-density 
regions into which the radiation streams. 

Although radiation does not prevent accretion entirely, it does heat the 
surrounding core once a protostellar object has formed in the center. The lu- 
minosity from accretion onto even a sub-solar mass protostar can heat the sur- 
rounding envelope enoug h to prevent fra gmentation during the collapse of a 
centrally condensed core (jKrumholzl 120061 ) . Adaptive mesh refinement simula- 
tions including radiative transfer using a temperature dependent opacity show 
that fragmentation is strongl y suppressed when compared to models with an 
isothermal equation of state (jKrumholz. Klein &: McKed 120071 ). At the time of 
publication, these models had only reached 10 M Q , but more recent conference 
reports have not yet given an upper mass li mit. Although this work calls int o 
question the competitive accretion scenario (jKrumholz. Mckee fe Kleinl l2005bl ). 
it remains to be seen whether the two scenarios converge to a common descrip- 
tion of collapse when beginning from a turbulent cloud rather than a centrally 
condensed core (see Clark or Bonnell in this volume). 

The good news here is that we seem to understand why feedback does not 
entirely prevent massive star formation. The bad news is that we still don't 
know if it does or does not ult imately determine the upper mass limit of around 
100-150 M (e.g., Figer 2005). Several alternate explanations can be imagined. 
For example, if fragmentation cannot be prevented for cores above the upper 
mass limit (perhaps because of insufficient central radiative heating), then the 
upper mass limit would be determined by the largest mass reservoirs available. 
Alternatively, disk fragmentation could cut off accretion around too massive 
stars. Another suggestion is that massive stars grow by direct collision, so that 
the stellar density is the limiting factor. 

If feedback does act to cut off accretion, which mechanism is dominant? Jets 



and w inds appear to evolve during the formation of massive stars. iBeuther Shepherd 
(2005) suggest an evolutionary sequence in which the outflow type evolves as the 
central star accretes additional mass. In this sequence, high mass protostellar 
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objects with central stars appearing as mid-to-early B stars have strongly colli- 
mated jets, early B to late O stars that have already formed hypercompact H n 
regions start to show a line-driven wind in addition to the jet, while early O 
stars that have accreted most or all of their mass and have ultracompact H n 
regions (UCHRs) lose the jet and retain only the wind. This sequence could 
be symptomatic of disk destruction and accretion cutoff in the final stages, 
a hy pothesis supported by the lack of observation of Q sta r s with collimated 
jets (IShepherd. Test. & Stark] 120031 : ISollins fc Meeeathl 120041 : lArce et al.1 120071 ) 



or disks ( 



Cesaroni et al.l 120071 ). Whether the winds themselves can destroy the 



disks remains an open question, however. 

Jets have been suggested bv lNakamura TA ( 20071 ) as a way of supporting 
cluster-forming cores. However, the relatively low-resolution numerical models 
they used do not yet resolve the question of whether the turbulence driven in the 
cluster can actually cut off accretion. Furthermore, their outflow-driven turbu- 
lence model predicts a prominent break in the kinetic energy power spectrum at 
a scale close to the outflow len gth. Such a break has not yet been observed in a 
molecular cloud. For example. lOssenkopf &: Mac Low! ( 2002) found no b r eak in 

(120071 ) 
and find 



Baneriee et al 



the Polaris Flare — admittedly not a cluster forming core 
examine the velocity distribution produced by a single Mach 5 jet 
very little supersonic material outside the head of the jet, suggesting that jets 
may have diffic ulty driving the observed supersonic turbulence in clouds. On 
the other hand, Matzner ( 20071 ) offers an analytic treatment that argues that 
j ets ca n indeed drive the observed turbulence, consistent with iNakamura &: Li 
(120071 1 . 

Ionizing radiation fr om th e growi ng central star offers another mechanism 
for cutting off accretion. Keto ( 2002b! ) made the important point that the for- 
mation of an H II region around a massive star does not cut off accretion until 
the Stromgren radius grows to beyond the Bondi-Parker radiu s wh e re the sound 
speed in the ionized gas exceeds the local escape velocity. iKetd (l2007f ) semi- 
analytically predicted the evolutionary sequence for star s forming in r o tating , 
collapsing, cores, using the similarity solution derived by iTerebev et al.l ((1984). 
A gravitationally confined, quasi-spherical, hypercompact H n region forms first, 
which then transforms to a bipolar UCHR as accretion increases the central 
star's mass and luminosity. As the ionizing luminosity increases further, only 
a remnant disk remains neutral, and finally that evaporates as well. The basic 
prediction of this model that hypercompact H n regions should sh ow ey i dence o f 
accreting ionized gas has been borne out by observations both bv lKetol (2002a); 
iKeto & Wool (120061 1 and by others (jBeltran et al.ll2006h . 

If ionizing radiation is to cut off accretion, the final stage of accretion will 
be through a photoevaporatin g disk. This configuration was first analytically 
described by iHollenbach et al. (1994) in the context of solar mass star forma- 
tion, while in the context of massive stars a n analytic treatment of the str ucture 
of the photoevaporating wind was done b ylLugo. Lizano. h Garayl ( 20041 ). The 
best numerical model to date was done by iRichling &: Yorkel (jl997l ) using a two- 
dimensional radiation gas dynamics code that included an explicit treatment of 
dust scattering. Dust scattering leads to a factor of 3-4 increase in the photo- 
evaporation mass loss rate, which they show to be p roportional to th e stella r 
luminosity <S*' 58 , close to the analytic prediction of Hollenbach et all ( 19941 ). 
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though flattening somewhat at photon luminosities 5* > 10 47 s _1 . They also 
explicitly included line-driven stellar winds in some models, demonstrating that 
they were collimated by the disk and photoevaporating wind. 

If UCHRs expand at the sound speed of ionized gas, c« ~ 10 km s , they 
should have lifetimes of roughly 10 4 yr. Less than 1% of an OB star's lifetime of 
a few megayears should therefore be spent within an UCHR, so the same fraction 
of OB stars should now lie within UCHRs. However. IWood & Churchwelll (|1989l ) 



surveyed UCHRs and found numbers in our Galaxy consistent with over 10% of 
OB stars being surrounded by them, or equivalently, lifetimes > 10 5 yr. 

A number of explanations have been proposed for this lifetime problem, 
including thermal pressure confinement in cloud cores, ram pressure confinement 
by infall or bow s hocks, champagn e flows, disk evaporation, and mass-loaded 
stellar winds (see Churchwelll Il999l ). Several of these explanations have basic 



problems that suggest they likely cannot explain the lifetime problem. 

Confinement by thermal pressure of t he surrounding molecular gas re- 
quires pressures of P/k y 10 8 -10 9 cm" 3 K ( De Pree. Rodriguez. Gossi[T995l : 



Garcia-Segura Franco! fl996l ). At typical molecular clou d temperatur es of 10 
100 K, this implies densities n > 10 6 cm" 3 . However, the Jeansl ( 19021 ) mass 



Mj = {lf) fe) = 61M °{w^) (wok) ' (1) 

where we have assumed a mean mass per particle fi = 3.87 x 10~ 24 cm -3 ap- 
propriate for fully molecular gas with one helium atom for every ten hydrogen 
nuclei. Therefore cores massive enough to form OB stars c ontain multiple Jeans 



mass es and are thus very likely to be freely collapsing (e.g. iMac Low fc Klessen 
20041 ). The free-fall time 



^=(32^) = (3-4 x 10 yr) \JJ^~^~3j • (2) 

Typical lifetimes of > 10 5 yr would thus require massive cores to last > 3tg at the 
hypothesized densities, rather than dynamically collapsing. Although these high 
pressures are indeed observed, the y are unlik e ly to occur in objects with lifetimes 
long enough to solve the problem. IXie et ail ( 19961 ) instead proposed a variation 



on this theme: confinement by turbulent rather than thermal pressure. However, 
turbulent motions decay quickly, with a characte ristic timescale of less than 
a fre e -fall time unde r molecular cloud conditions ( Stone. Ostriker. &: Gammid 



1998: IMac Lowl H999). Turbulent pressure would thus have to be continuously 
replenished to maintain confinement for multiple free-fall times, which would be 
difficult at such high densities and small scales. 

Another option is ram pressure confinement of UCHRs by infall of sur- 
rounding gas. However, this is unstable for two different physical reasons. First, 
the density and photon flux will follow different power laws in a gravitationally 
infalling region i onized from within, so they can neve r balance each other in sta- 
ble equilibrium ( Yorkd[l986l : Hollenbach et al. 19941 ). Either the ionized region 



will expand, or the infall will smother the ionizing source. Second, the situation 
is Rayleigh- Taylor unstable, as this option requires the rarefied, ionized gas of 
the UCHR to support the infalling, dense gas. 
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Bow shock models ( Van Buren et al.1ll990l ; lMac Low et al.lll99ll : lArthur Sz Hoare 



a vc- 



2006) require high values of ram pressure P ram oc nvj- At n = 10 5 cm 
locity of ~ 10 km s _1 is required (jVan Buren et aL 19901 ). A star moving at such 
a high velocity in a straight line would travel a parsec over the supposed UCHR 
lifetime of 10 5 yr, requiring a uniform-density region of mass > 5 x 10 3 M for 
confinement. As collapse would occur on the same timescale, more mass would 
actually be required to solve the lifetime problem. 

Expansion of an H n region in a density g radient can drive su personic 
champagne flows down steep enough gradients ( Tenorio-Taglei Il979l ). Two- 
dimensional models first studied expansion a cross sharp density discontinuities 
( Bodenheimer. Tenorio-Tagle. Yorkd 19791 ). but then examined other config- 
urations such as freely-collapsing cloud cores (lYorke et al.lll982D. clouds with 
powe r-law density grad ients in spherica l ( Franc o. Tenorio-Tagle. fc Bodenheimer 
1990), and cylindrical ( Garcia-Segura & Francol ll996) configurations, and expo- 



nential density gradients ( Arthur Hoaref 2006 ) . Stellar wind combined with a 
champ agne flow down a continuous gradient was also modeled bv lArthur Sz Hoard 
( 20061 ) . This model does fit well the observed velocity structure in some cometary 
UCHRs such as G29.96— 0.02. However, these models face the same timescale 
problem as thermal pressure confinement models: regions dense enough to 
explain the observations are gravitationally unstable, and collapse on short 
timescales. 

Another class of models re lies on mass-loaded stellar winds to reproduce the 

observed properties of UCHRs ( Dyson. Williams. Redmanlll995l ; iRedman. Williams. Sz Dyson 
1996j; IWilliams. Dyson. Sz Redman 1996 ; Lizano et al. 19961 V In these mod- 



els, an expanding stellar wind entrains a distribution of small, self-gravitating, 
pressure-confined clumps that take substantial time to evaporate. These models 
can repr oduce many of the basic features of the observations incl uding some line 



profiles (|Dvson et al 



and bipolar shapes ( 



1995 ), core-halo, shell ( Redman et al.l 19961 ) , and cometary 



Redman. Williams. Sz Dyson 19981 ). but at the cost of re- 



quiring an arbitrary distribution of pre-existing clumps that cannot be self- 

consistently predicted. 

Finally, iMac Low etaD (|2007l ) examines the expansion of H n regions into 
turbulent, self-gravitating gas. At densities high enough for massive stars to 
form, the expandi ng shell dr i ven by newly ionized gas qu ickly becomes gravita- 
tionally unstable (Voit 19881 ; Mac Low Sz Norman 19931 ). collapsing even more 
promptly than the surrounding gas. These regions of secondary collapse in the 
shell may be externally ionized to form objects with the properties of UCHRs. 
As the shell expands to larger sizes, new regions can form, extending the lifetime 
during which UCHRs remain visible well beyond the expansion time of the orig- 
inal H II region. Some combination of secondary collapse and disk evaporation 
appears to me most likely to ultimately be the solution to the UCHR lifetime 
problem. 



2. Large Scale Feedback 

The radiation, winds, and supernovae from massive stars also influence their 
galactic environment. Two questions of current interest are: 1) Can compres- 
sion driven by feedback trigger further star formation with an efficiency above 
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unity, allowing self-propagating star formation? 2) Can any of these forms of 
feedback support molecular clouds for many dynamical times, reducing their star 
formation efficiency per unit free-fall time, or, alternatively, giving a negative 
triggering efficiency? 

Stellar feedback certainly produces compressive motions. The suggestion 
that expanding H n regions might compress s urrounding gas and t rigge r fur- 
ther massive star forma t ion w as first made by lElmegreen Sz Lada ( 19771 ) and 
Elmegreen Sz Elmegreen ( 19781 ). They assumed that turbule nt velocities in the 



shell w oul d be of the sam e order as the expansion velocity. lOstriker Sz Cowiel 
( 198ll ) and lVishniad ( 19831 ) assumed, on the other hand, that the turbulent veloc- 
ities in the shell woul d onl y be tr ansonic. This latter assu mption was supported 
by an alytic work by IVoitl ( 19881 ) and numerical work by iMac Low Sz Norman 
( 19931 ) . The operation of this mechanism in three- dimensional simul ations of 
expansion into a turbulent medium is describ ed by Dale et al. ( 20051 ) using a 
smoothed particle hy drodynamics c ode and by IMac Low et ah ( 20071 ) using the 



ZEUS-MP grid code ( Norman 20001 ). These models show that the dumpiness of 
the turbulent medium allows the ionization front to blow out of the collapsing 
core, while the denser regions resist ionization and continue to collapse. Neither 
m odel finds evidenc e for efficient triggering — indeed in the more global model 
of iDale et all (|2005T ) there is some evidence for negativ e feedback — fewer stars 
formin g in the presence of the H II region than otherwise. IDale. Clark. Sz Bonnell 
(2007) extended this result by examining an externally ionized collapsing core, 
and found rather inefficient triggering in that case as well. These low or negative 
triggering efficiencies occur when the acceleration and stirring of the gas by a 
shock wave expanding into a clumpy medium dominates over direct compression. 

Supernovae compress gas on even larger scales than H n reg i ons, so they 
are also often invoked as triggering agents. Joung Sz Mac Low ( 20061 ) simu- 
lated surjenrova^riveriturbulence using the Flash adaptive mesh refinement 
code iFrvxell et all (|2000l ) on a grid covering 500 pc on a side in the plane and 
5 kpc up and down vertically, with 2 pc resolution within 200 pc of the plane. 
When the galactic fountain flow had reached a steady state, the density, pres- 
sure, and velocity dispersion of the gas in two-zone cubes was compared to the 
turbulent Jeans criterion for collapse, to measure the expected star formation 
rate in the flow. If a generous local star formation efficiency of 30% is assumed 
for the unstable gas, they find that the resulting star formation rate is a full 
order of magnitude below the value required to generate the input supernova 
rate. 

These results on triggering by ionizing radiation and supernova explosions 
taken together suggest that triggering does indeed occur, but that it is only a 
10% effect rather than the primary mechanism controlling star formation (al- 
though it may be more important than this for massive s tar formation i n par - 



ticular). A recent observational test supports this result. Mizuno et al. ( 20071 ) 



used the NANTEN telescope to survey dense molecular gas interacting with 23 
southern H n regions within 4 kpc of the Sun. The detected clouds were divided 
into regions adjacent or not adjacent to the H n regions. The adjacent regions 
indeed contained more protostellar objects, and ones with significantly higher 
far- infrared luminosities luminous protostellar objects. However, when the effect 
is integrated over the galaxy, it appears to be only a 10-30% effect. 
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If triggering is not that effective, is feedback good at the converse — supporting 
molecular clouds against collapse? The question of how long-lived molecular 
clouds remains very controversial, with answers ranging over an order of magni- 
tude from one to as long as te n free-fall times. Long time scales a re argued for 
on the basis of either magnetic (iMouschovias. Tassis. fe Kunz 2006) or turbulent 
(jKrumholz. Matzner. McKed I2006T) support, while short time scales depend 



on sweeping by shock compressio n (Ballcsteros-Parcdes, J. Hartmann . L.. Vazquez-Semadeni. E. 1999. Ap 



19991 ; Hartmann. Ballesteros-Paredes. Berginll200ll; Padoan et al.ll2007h or di 
rect g ravitational instability and collapse (ILi. Mac Low, fc Klessenll2005l ; lElmegreen 



2007) 



iMatzner ( 20021 ) argued that turbulent support of molecular clouds can 
be maintained by H n region expansion within the clouds, for clouds with a 
mass greater than 3.7 x 10 5 Mq. His analysis suggests that such large clouds 
reach a balance between energy inj ection and decay in a turbulent cascade. 
Krumholz. Matzner. McKed ( 20061 ) used this analysis and made the additional 



assumption that giant molecular clouds are spherical objects with homologous 
power-law distributions of density, energy injection, and oth er relevant quanti- 
ties. Using a star formatio n law dKrumholz &; McKed I2005T) derived from the 
turbu lent collapse models of lVazquez-Semadeni. Ballesteros-Paredes. Sz Klessen 
(2003), they find that large molecular clou ds remain quasi - stable for as long as 
20-40 Myr. Similar work was performed bv lHuff &: Stahler ( 20071 ) in the context 
of the Orion Nebula Clus t er. 

However, Elmegreen ( 20071 ) reaches opposite conclusions working from al- 
most the same assumptions. He argues that giant molecular clouds never have 
spherical, homologous shapes, but rather are filamentary objects containing 
quasi-spherical cluster-forming cores of order 10 M . In his picture, these 
dense cores are dispersed by H n region expansion within 2-4 Myr, but mag- 
netically supported envelopes that are not participating in gravitational collapse 
may indeed survive for much longer. 

The observational evidence has not yet clarified the picture either. A 
study by Fukui's group of the associations between you ng clusters, H n re gions, 
and molecular gas observed by NANTEN, reported in Blitz et ail (2006), sug- 
gests that GMC s rema in without visible H n regions for about 7 Myr — though 
iGruendl Chul ( 20071 ) used Spitzer to show that they do have young stellar 
objects. About another 14 Myr passes before visible clusters appear, and then 
6 Myr remain before the molecular gas is cleared away, for a total of as much as 
27 Myr. They also note that the average mass of the observed objects increases 
with the stage of life. They interpret this as suggesting continuing accretion 
onto the clouds over their lifetimes. 

On the other hand, iTamburro et al.l ( 20071 ) have used g lobal imaging of 



nearb y galaxies from the Spitzer IR Nearb y Galaxies Survey (jKennicutt et al 



l2003h and The H I Nearby Galaxy Survey (jWalter et al.1 120071 ) to measure the 
average delay between the peak of H I emission and the peak of star formation 
as measured by 24 fim emission from warm dust in regions of massive star 
formation. They find values of 2-4 Myr for the delay, with no values above 
6 Myr in their sample of over 20 galaxies. 

While the peak of 24 /xm emission does not yield the total lifetime of the 
molecular clouds, it does measure the period of strongest star formation. The 
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results of iTamburro et all ( 20071 ) suggest that star formation begins promptly 



and strongly within a few million years of the on set of gravit a tiona l instability as 
traced by H I emission. The timescales cited bv lBlitz et al.l ((2006) make no ref- 
erence to the strength of star formation during any particular period, though, so 
perhaps the two observations are consistent with the theoretical picture of cloud 
support by H n region expansion: if clouds collapse promptly, forming many 
young stars, the result could be well supported clouds with low ongoing star 
formation efficiency. Investigation of this scenario looks like a fruitful direction 
for further research. 
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